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o": 

C> This paper derives integral representations for the Black-Scholes 

price of arithmetic-average Asian options. Their proof is by 
Laplace inverting the Laplace transform of [GY] using complex 
analytic methods. The analysis ultimately rests on the gamma 
function which in this sense is at the base of Asian options. The 
^5 ■ results of [GY] are corrected and their validitity is extended. 

o\ : 

i — 1| Introduction: The aim of this paper is to derive exact closed-form valuation formulas 
for European-style arithmetic-average Asian options in a Black-Scholes framework. These 
options are a particular class of path-dependent options with the arithmetic average of 
■ the underlying security as the relevant variable. 

This paper has its origin in that developement in the analysis of Asian options initiated 
by [Y]. Yor's valuation formula gives clear evidence that pricing Asian options is a problem 
of some intrinsic difficulty indeed for which no, in the strict sense, simple solution should be 
expected. Instead, one should, as a first step, ask for structurally clear solutions, and only 
then, as a second step, consider in particular questions of a more computational nature. 

in ■ 

The main valuation result of this paper is described in paragraph four and identifies the 
Black-Scholes price of an Asian option as an integral over, roughly speaking, the product of 
two well-studied higher transcendental functions both given as integrals and built up using 
so-called Hermite functions. These last functions generalize the familiar complementary 
error function and come from boundary-value problems in potential theory for domains 
whose surface is an infinite parabolic cylinder. It seems to be the first time that such 
parabolic cylinder functions are explicitly identified to characterize solutions to problems 

L" ■ of Finance. Their relevant properties are reviewed in paragraph zero. 

^ • 

Both valuation formulas involve three integrations. In a sense, this measures the difficulty 
of the problem. For getting some intuition, first think of the error integral as an example of 
how a single integration yields a higher transcendental function. In the best case, a triple 
integral separates as a product of three such functions. In the worst case, no separation 
is possible and all one can say is that one has to integrate a function in three variables 
three times in succession. With its integrand a function f(x,y,w) that cannot be further 
separated, Yor's formula belongs to this last category. In our formula complete separation 
is not achieved either. However, with its integrand a product g(x,y)h(y,w) it is a single 
integral of a function given as the product of two single integrals. 

The method used in this paper is the Laplace transform Ansatz to option pricing of [GY] 
and its main result is proved by Laplace inverting the Laplace transform they have derived. 

Transform methods like the Laplace transform are classical concepts in analysis with ap- 
plications ranging from number theory to boundary-value problems for partial differential 



1 



2 



Michael Schroder 



equations. In particular the Laplace transform approach used in this paper generally 
consists of two stages. As a first step, a suitably nice function on the positive real line is 
transformed into a complex-valued function. This Laplace transform exists on a half-plane 
sufficiently deep within the right complex half-plane and defines there a holomorphic, i.e., 
complex analytic, function. The methods in this part of the transform analysis have real 
analysis origin. As a second step, the Laplace transform has to be inverted to give the 
desired function on the positive real line. The natural methods in this step are complex 
analytic. Some relevant concepts are recalled in paragraph zero of this paper 

The computation of the Laplace transform already is quite a problem in general. Indeed, 
with the functions to be transformed generally not known one expressed aim for trying to 
calculate their Laplace transform actually is to get an explicit expression at all. However, 
for an actual calculation of a Laplace transform in such a situation additional insights 
into the mathematical structure are necessary. Thus, the Laplace transform in [GY] is 
computed introducing a stochastic clock and using general results on Bessel processes. 

It should be said that, unfortunately, there are serious problems of a fundamental nature 
with this result. More precisely, the conceptual error I found in September 1999 kills this 
approach to the valuation of Asian options as it stands. Luckily, Peter Carr eventually 
succeeded in convincing me that the Geman-Yor result can nevertheless be used in a very 
astute way for its original purpose. All this is discussed in paragraph five of this paper. 

Mathematically, the inversion of a non-trivial Laplace transform belongs to the nastier 
problems in analysis. Difficulties are encountered here as a rule, and more than often one 
has to be satisfied with only knowing the Laplace tranform. It thus can be regarded as as 
the main mathematical contribution of the paper that it indeed provides such an inversion 
for the above Laplace transform of the price of Asian options. The details of the proof 
are described in paragraphs six to twelve of this paper. Its methods are complex analytic. 
It uses the complex inversion formula for the Laplace transform and is ultimately based 
on the classical Hankel formulas for the gamma function. These methods are such that 
the results can be extended from the classical Black-Scholes framework essential for the 
validity of [GY] to one with no restrictions on the drift coefficients. 

It came as a surprise to the author that one has such a structural integral representation 
for the value of Asian options as derived in this paper. He is not sure if this should be 
attributed to the qualities intrinsic to Asian options or to the fact that, as described above, 
one takes a somewhat indirect approach to their valuation. 

The series representations and asymptotic expansions of our formula appropriate and ad- 
equate for computation are discussed in [Sch]. 
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0. Analytic preliminaries: This paragraph collects some relevant facts about the 
Laplace transform from [B], [D], and about Hermite functions from [L, §§10.2ff]. 

Laplace transform: The class of functions considered in the sequel is that of functions 
of exponential type, i.e., of continuous, real-valued functions / on the non-negative real 
line such that there is a real number a for which exp (at)f(t) is bounded for any t > 0. The 
Laplace transform is the linear operator L that associates to any function / of exponential 
type the complex-valued function L(/), holomorphic on a suitable complex half-plane, 
which for any complex number z with sufficiently big real part is explicitly given by: 



POO 

L(/)(z)= / e- zt f(t)alt. 
Jo 



L 



This operator is an injection. Its inverse, the inverse Laplace transform L -1 , is expressed 
as a contour integral by the complex inversion formula of Riemann. Applying to any 
function H analytic on half-planes {Re (z) >zq} with zq any sufficiently big positive real 
number, it asserts: 

-\H){t) = ^- e zt -H(z)dz, 

for any positive real number t if H satisfies a growth condition at infinity such that the 
above integral exists. For any function / of exponential type one so has in particular 
(L- 1 oL)(/) = /. 

For getting an idea of the proof of the complex inversion formula start with any function / 
holomorphic on a half-plane {z|Re(z) > a}. For obtaining a suitable integral representa- 
tion of /, invoke the Cauchy integral formula in the following way. Fix z with Re (z) > a 
and choose a line C(z ) = {z\Re (z) = z } such that z is to the right of it, i.e., z < Re (z), 
and such that it is contained in the half-plane where / is holomorphic, i.e. a < zq. Make 
this line into a path of integration by moving upwards from z$ — zoo to zq + zoo. Shift pe- 
spective to the Riemann sphere by adding the point infinity to the complex plane. In this 
picture, C(zq) is a closed path that circles clockwise, i.e., in the mathematically negative 
direction, around z. Formally apply the Cauchy integral formula, expressing 2%i ■ f(z) as 
minus the integral of (w — z)~ 1 f(w) over C(zq). Since — (w — z)~ l is the Laplace transform 
at z of the map on the positive real line given by exp (wh) for any h > 0, it so follows: 



f(z) = — [ [ e {w - z ^ h f{w)dhdw 

2m JC(zn) JO 



'C(zo) JO 

On the Riemann sphere the contour C(zq) can be seen as the limit of closed contours 
contained in the complex plane. The above integral thus makes sense as such a limit if 
/ is supposed rapidly decreasing, i.e., the absolute value of / exponentially decreases to 
zero with the absolute value \w\ of w going to infinity. Then apply Fubini's theorem to 
interchange the integrals in the above formula. The inverse Laplace transform of / is then 
identified as the function on the positive real line given by 



1 

2ni 



e wh f(w) dw 

C(z ) 



for any h > 0, which is the complex inversion formula. 
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Hermite functions: The Hermite function H v of degree any complex number v is the 
complex-valued function on the complex plane given as solution to the differential equation 
u" — 2zu' + 2vu = in the complex variable z. It is holomorphic on the complex plane as 
a function of both its variable z and its degree v. Using the differential equation, for v any 
non-negative integer, YL V is the z/-th Hermite polynomial and belongs to a standard class 
of orthogonal functions. If the real part of the degree v is negative the Hermite function 
H v has the following integral representation: 

H U { Z ) = — — / e -v 2 -2zu u -(„+l) du ; 

1 \— v ) Jo 

for any complex number z and with r(— v) the value of the gamma function at —v. Her- 
mite functions thus also generalize the complementary error function Erfc . Recalling that 
Erfc is for any complex number z given by any of the following expressions: 

ry POO ry pCO 

Erfc (z) = N(-z) = — / e- u2 du = e" 2 ' -= / e~^- 2xu du , 

V 71 " Jz V 71 " Jo 

notice (2/y/n) ■ H-i(z) = exp (^ 2 )Erfc (z). More generally, one has for any non-negative 
integer n the identity: 

, 0F (-1)" d n ( z 2 x 

H_ (n+1 ) (z) = ^— — — e z Erfc (z) 

{n+ n ' 2 nl ■ 2 n dz n V v ') 

for any complex number z. Hermite functions of degree any complex number v satisfy the 
recurrence relations: 

H' v (z) = 2vH v - 1 (z), 

H v+1 (z) - 2z ■ H v (z) + 2v ■ H u -!(z) = 0. 

Using the differential equation of the Hermite functions the second of these is implied by 
the first which can be seen using the following absolutely convergent series expansion: 

valid for any complex number z and uniformly convergent on compact sets. Like the ex- 
ponential series this series has excellent convergence properties for complex numbers z of 
smaller absolute value. For arguments z of larger absolute value one uses asymptotic ex- 
pansions. If the real part of v is negative the Hermite function of degree v has the following 
asymptotic expansion in the right half-plane: 

for any complex number z with positive real part and any non-negative integer n. Herein 
(— v) 2 k is the Pochhammer symbol recalled to be given by (— u) = 1 and (— v)2k = 
r(— v + 2fe)/r(— v) = (—v)(—v + 1) • • ■ (— v + 2k — 1) for k positive. For estimating the 
error term r n +i, fix any positive real number 5 less than or equal to tv/2 and choose any 
complex number z in the wedge contained in the right half-plane that is enclosed between 
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the two rays emanating from the origin with angles ±(7r/2 — 5) respectively. Equivalently, 
the absolute value of the argument of z is at most n/2 — 6. Then the absolute value of the 
error term r n+ i(z) satisfies the following estimate: 

I ( )\< 1 r(2(n + l)-Re(z/)) 1 

lTn+l[Z)l ~ ( sin ( 5 )) 2 ^ +1 )- Re H ' (n+l)\\T(-v)\ (2|z|)2(«+l)-ReH • 

The differential equation for the Hermite function of degree any complex number v can be 
transformed into the differential equation u" + {y + 1/2 — z 2 /4)u = whose solution D v is 
thus connected with the Hermite function H v by D v (z) = 2~ v l 2 exp (— z 2 /4)H v {2~ l / 2 z), 
for any complex number z. Similarly, Hermite functions are related to the second Rum- 
mer confluent hypergeometric function \I> by H u {z) = 2-\I/(— v/2, 1/2, z 2 ), for any complex 
number z in the right half-plane. 



Part I Statement of results 

1. Black— Scholes framework for valuating contingent claims: For our analysis 
we place ourselves in the Black-Scholes framework and use the risk-neutral approach to 
valuating contingent claims as described in [DSM, Chapters 17, 22]. 

In this set-up one has two securities. First there is a riskless security, a bond, that has the 
continuously compounding positive interest rate r. Then there is a risky security whose 
price process S is modelled as follows. Start with a complete probability space equipped 
with the standard filtration of a standard Brownian motion on it that has the time set 
[0, oo). On this filtered space one has the risk neutral measure Q, a probability measure 
equivalent to the given one, and a standard Q-Brownian motion B such that S is the 
strong solution of the following stochastic differential equation: 

dS t = w ■ S t dt + a ■ S t dB t , te[0,oo). 

Herein the positive constant a is the volatility of S, whereas the specific form of the con- 
stant w depends on the security modelled. For instance it is the interest rate rate if S is 
a non-dividend paying stock. 

With the existence of Q one has the arbitrage-pricing principle: let Y be any European- 
style contingent claim on the above filtered probability space written at time to an d paying 
Yt at time T. At any time t between to and T its time-t price C t is then: 

C t = e- r <- T -*E?[Y T ] 

with the expectation conditional on the information available at time t and taken with 
respect to Q. 

2. The notion of arithmetic— average Asian options: Fix any time to an d consider 
in the Black-Scholes framework reviewed above the process J given for any time t by: 

J(t) = f S T dr. 

J to 



Michael Schroder 



The European style arithmetic-average Asian option written at time to, with maturity the 
time T, and fixed strike price K is then the contingent claim on [to, T] paying 



(J{T) 



K 



+ 



:= max < 0, 



J(T) 



K 



\T-t ") ~ T-to 

at time T. Recall that points in time are taken to be non-negative real numbers. 

3. The valuation problem for arithmetic— average Asian options: Using the 
arbitrage-pricing principle, the price Ct,r{K) at any time t between to and T of the 
arithmetic-average Asian option introduced in paragraph two is given by: 



C t , T {K) = e-« T -*E? 



J{T) 



K 



T-to 

From [GY, §3.2] one has the following normalization of this value process: 

Lemma: For any time t between to and T , one has: 

e-r(T-t) 4 o 



where 



C t , T (K) 



C {u) (h, q) = E Q 



T-to 



max jo, jf e 2 ^ + ^dr-qj . 



To explain the notation, the normalized parameters herein are as follows: 

2vj 



v 



- 1, 



h=-(T-t), 



q =^(K. { T-to)-j\ d r). 



To interpret these quantities, C^(h, q) is the normalized price for the Asian option. It de- 
pends on the normalized adjusted interest rate v, which for w positive is bigger than minus 
one, on the normalized time to maturity h, which is non-negative, and on the normalized 
strike price q. 



To prove the Lemma break up the integral J(T) at the point t and bring St in front of 
the expectation. Modify that part of the integral unknown given the information at time 
t as follows. Change variables to obtain an integral from time to time T—t. The strong 
Markov property of Brownian motion identifies the exponent B(t + r) — B(t) of the new 
integrand as a standard Brownian motion that is independent of the information at time 
t. Herein change time by dividing through a 2 / 4. Using the scaling property of Brownian 
motions, this time-changed process is again a standard Brownian motion if multiplied by 
a/2 which completes the proof. 
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4. Statement of results: In the setting of paragraphs one to three consider any time 
t between the times to and T. Using the decomposition 

-r(T-t) 40 

of §3 Lemma, the valuation of the arithmetic-average Asian option is reduced to comput- 
ing its normalized price C^ u \h, q). This price depends on the sign of the normalized strike 
price q. This parameter is recalled to be given by: 

a 2 ! / r T-(4/a 2 )-h 



[K-(T-t )- / S u du). 

i/a 2 )-h V Jt () J 



T-{A/ai)-h \ J t 

as a function of time t which enters via the normalized time to maturity h = (a 2 / 4) ■ (T—t). 
Generically q will be positive. It is non-positive if and only if at time t it is already known 
that the option will be in the money at its maturity T. 

The main result: The main result of this paper is a closed form solution for C^\h, q) 
if q is positive. It expresses this function as the sum of integral representations. These 
are obtained by integrating the product of Hermite functions, discussed in paragraph zero, 
with functions derived from weighted complementary error functions. Proved later in 
paragraphs nine to twelve, the precise result is as follows: 

Theorem: If q is positive, the normalized price C("\h,q) of the Asian option is given 
by the following five-term sum: 

I>+4) (2q)^ / x 

g 2tv(u+1) ' ~ ' [ Lyh yP^+ 2 ~+~ U hyp,-(i/+2) - ^hyp,;/ ~ ^hyp,-^ J ■ 



e 



To explain the functions, recall from paragraph zero the Hermite functions of degree 
any complex number fx. The trigonometric term C tr i g is then: 

Ctrig = c • j H_( u+4) (^- C °^) ^j cos(v -6)d0, 

recalling v = 2zu/a 2 — l and abbreviating: 

c = c(v, q) = I>+4) • • . 

u+1 n-e^ 

For v bigger than minus one C tr i g has the additional representation: 



e 2h(u+i)_ 1 „oo /coshfxU 



— vx , 

e ax . 
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The hyperbolic terms Ch yp ,& with b equal to ±u or ±(i/+2) are: 

= ^ ' e^ + ^ jf° H_ {u+A) (^^) • • E b (h)(y) dy . 
Herein E b {h) is in terms of real-valued functions for any real number y given by: 



i du 



It is identified in paragraph ten as imaginary part of a certain weighted complementary 
error function at complex arguments. Actually, this last interpretation seems to be best 
suited for a numerical analysis. 

The case q < 0: It is a characteristic feature of Asian options that their prices become 
rather rigid some time before maturity. The present case q < where at time t it is known 
that the option will be in the money at its maturity makes this precise. The rigidity of 
the price is reflected in the formula of [GY, §3.4]. For any h > one has: 

e 2(u+l)h _ 1 

The proof reduces to computing the Q-expectation of J Q h exp (2(£? T +z/r)) dr. For an argu- 
ment independent of [GY, §3.4], applying Fubini's theorem this expectation is the integral 
from zero to h with respect to the variable r of the product exp (2z/r) times the expectation 
of exp(2S T ). This last expectation is the integral over the real line with respect to the 
variable x of the product exp (2x) times the density (27rr) -1 / 2 exp (— (2t)~ 1 x 2 ). Chang- 
ing variables y = (2r) _1 / 2 (x — 2r) it is seen to be equal to exp (2r). On substitution of 
this result one so is reduced to calculate the integral from zero to h with respect to r of 
exp (2(z/ + 1)t). Distinguishing the cases where v + 1 is zero and non-zero, this is seen to 
have the above value, completing the proof. 

Remark: The author is very grateful to a number of people, including a referee of this 
paper, for having drawn his attention to Yor's formula [Y, (6.e), p. 528]. His integral 
representation is a triple integral with a trigonometric-hyperbolic type integrand: 

CM(h,q) = j==- \ / e 2 -(--q) -^ xy (h)dydx, 

7rV2/i Jo Jo K v 1 

where for any positive real number a, the function ifj a is given for any h > by: 

/ /M r ~ih -«cosh(«,) . /7T \ 

iPa{h) = J e e smh(w) ■ sm y—w J dw . 

In contrast, our formula separates into a trigonometric and a hyperbolic part. It is given as 
a sum of single integrals whose integrands have a structural interpretation as products of 
two functions. It moreover identifies the higher transcendental functions occuring as factors 
in these products, and shows how they are given by or built up from Hermite functions. On 
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a technical level, these differences can be regarded as consequences of the different mathe- 
matical approaches for proving the valuation formula. Which of the above representations 
is better suited for which purpose in which situation remains to be analyzed. 



Part II Proof of the valuation formula 



5. Computing the price of the Asian option using Laplace transforms: The 

basic idea for the Laplace transform approach to valuating Asian options is not to consider 
a single above Asian option. Instead, take the exercise time T of the option as variable 
and consider today at time t all above Asian options whose exercise times T range from 
t to points in time very far in the future. Then, one wants to average the values of these 
options. This, however, is not possible in a naive way. One has to suitably weight these 
values and only then compute their average. This is effected by the Laplace transform. 

In the generic case with q positive a candidate for the Laplace transform of the factor 
CM( ,q) of the price for the Asian option was computed in [GY, pp.361ff]. Unfortu- 
nately, I found a conceptional error in their computations. Luckily however, Peter Carr 
eventually suceeded in making me see that this error is not as fatal as I originally believed. 
Indeed, in his interpretation the results of [GY] are sufficient for valuating Asian options 
as follows. 



Proposition: Suppose v = 2w / a 2 — 1 > — 1. If the positive real number h is such that 
q(h) = kh+q* is positive, the normalized price of the Asian option at h is given by: 

C<"> (h,q{h)) = L 1 (f gy {q(h),z)) (h) 

as the Laplace inverse of Fgy{q{^)i z ) °>t h. 



Leaving the function q unexplained for a moment, for any positive real number a define: 

D u (a,z) 



F GY (a,z) = 



z- (z-2(z/+l)) 



for any complex number z with positive real part bigger than 2(u+l), where on choosing 
the principal branch of the logarithm: 

D v (a,z)=?-p- [ e-^-x^-I^^^) dx . 

a Jq \a/ 

Here 1^ is, for any complex number //, the modified Bessel function of order fi more fully 
discussed in paragraph six below or in [L, Chapter 5], for instance. 

What regards the function q, the first general observation is that T does not enter directly 
in the analysis but indirectly via the normalized time to maturity h = (a 2 /A) ■ (T — t) of 
paragraph three. Letting T vary from t to infinity, h thus varies from zero to infinity. Take 
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h as new variable. The crucial difficulty now is that q of paragraph three also depends on 
h as a function of T. Indeed, abbreviating k = K / S t , on chasing definitions one gets: 

q(h) = k-h + q* where q* = ^— ( K ■ (t-t ) - [ S T dr 

4b t V Jt 

Herein q* is independent of T and thus of h. The upshot is that q is affine linear in h 
and non-constant in particular. This is a nasty situation since the computations of the 
normalized value of the Asian option in [GY, pp.361ff] erroneously take q as a constant. 
Indeed, they compute for any positive real number a and any complex number z with 
sufficiently big positive real part the Laplace transform: 

/•oo 

F GY (a,z):= / e~ zx -E[A^-a\dx, 
Jo 

setting: 

for any non-negative real number x, 

The proof of the Proposition then proceeds in two steps. Reduce to the above Geman-Yor 
Laplace transforms. Then sketch the essential steps of their computation. 

The basic idea of Peter Carr for reducing the Proposition to the Laplace transforms of 
[GY] was to fix any positive h and introduce an auxiliary function that coincides with the 
normalized value of the Asian option at h. More precisely, recalling §3 Lemma one has to 
consider the function ip given for any positive real number x by: 

tj}{x) = C (u) (x,kx + q*) = E[AP -kx-q*] . 

Deform this function into the function given for any positive real number x by: 

4>(x) = E[AP - kh-q*} . 

Then ip and are identical at h by construction. Granted it exists, the Laplace transform 
of (p is equal to FGY(q(h),z) recalling q(h) = kh+q*. Its Laplace inverse is the function (p 
by definition, and so equals ip when evaluated at h. One so has tautologically: 

CM (h, q(h)) = iP(h) = <j>{h) = L" 1 (F GY (q(h), z)) (h) , 

and is reduced to have the above Geman-Yor Laplace transforms Fgy{o,, z). 

The basic idea for computing these Laplace transforms in [GY, pp.361ff] is to introduce a 
suitable stochastic clock for the process J of paragraph two. This is made possible by the 
Bessel factorization result, attributed to Williams in [Y, §1.5] and to Lamperti in [RY, 
XI 1.28, p. 432], that the exponential of any Brownian motion with drift is a time-changed 
Bessel process in the following way: 



f X 2-{vu+B u ) 

= I e au , 

Jo 

to adopt Yor's notation. 
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Here R^ v ' is the Bessel process on [0, oo) with index z/, starting at with the value 1. At 
the level of a definition, the square of this process is a continuous diffusion process p with 
values in the non-negative real numbers satisfying: 

d Pt = 2(is+l)dt + 2^p- t dB t , p = l. 

To be able to compare times, let t UjU for any non-negative real number u denote the least 
upper bound of all non-negative real numbers s such that Jq exp (2(£? T + z/t)) <ir = u. 
Notice that this equals ^ (R^\s))~ 2 ds. 

The computation of Fqy (a, z) is then by deftly using the double role of as both the 
stochastic clock in its Bessel factorization and the actual underlying of the option. Indeed, 
fix any positive real number x, and consider the process at x on the set of all events 
where the passage time r v ^ a takes values less than or equal to x. Break the integral defing 
A^ v \x) at r VyCL . The first summand then is at time T Vja and so is equal to a. In 
the second summand one wants to restart the Brownian motion in the exponent of the 
integrand at r^ a . Thus shift the variable of integration accordingly. The second integral 
then is the product of exp (2 • (S(r !/ja )+z/ ■ T Uja )) times A^ u > at x—r v ^ ai by abuse of language 
after having applied Strong Markov. This last process is such that it is independent of 
the information at time r Uja . Recalling the role of r Uja , the first factor is the square of the 
Bessel process at time a. Now taking the expectation conditional on the information 
at r va , one thus gets: 



E Q 



(#)(i)-o) 



+ 



(R^\a)) 2 -E^ 



A i») 



x 



The Q-expectation of A^\w) = J W exp (2(B u +vu)) du is (exp (2(v+l)w) - l)/(2(i/+l)), 
as shown in §4 or using the more general results of [Y92, pp.69ff]. It thus follows: 



E Q {A^-ay = E® {R { "\a)) 



2(i/+l)[a:— tv, ] + 

2 e — 1 



2(i/ +1) 



To avoid a direct calculation of this expectation adopt the common strategy to first cal- 
culate its Laplace transform at arguments z whose real parts are positive and sufficiently 
big. The hope is to thus arrive at a simpler situation from which one is able to identify 
the original function itself. There is a technical point in that one wants to interchange 
the expectation with the Laplace transform integral. If z is real it is at this stage possibly 
best to follow Yor's proposal for justifying this. Indeed with the integrand of the double 
integral in question positive and measurable, apply Tonelli's theorem now but justify only 
in a later step that any of the resulting integrals is finite. The case of a general argument 
z is reduced to this case considering the absolute value of the integrand, and the result is: 



F GY {a,z) 



(R^(a)) 2 } 



z- (z-2(u+l)) 



One has enough information about Bessel processes for computing the Q-expectation of 
the right hand side. Using [Y80, Theoreme 4.7, p. 80], the expectation of exp (—z ■ T Vja ) 
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conditional upon BS U > (a) being the positive real number w is given by the following quotient 
of modified Bessel functions: 



L J I„ \a / 



This result ultimately depends on the following standard fact about Bessel processes [Y80, 
(4.3), p. 78]. If v is bigger than minus one, the density p Uja (l,w) of the Bessel semigroup 
with index v and starting point 1 at time a is given by: 



It is for this result that the hypothesis v bigger than minus one of the Proposition is 
crucially required. The numerator of the above Laplace transform now is obtained by 
integrating the product of these last two expressions and w square with respect to w over 
the positive real line. With the resulting integral finite for any complex number z with 
real part bigger than 2(z/+l), this completes the proof of the Proposition. 

6. Preliminaries on integral representations of gamma and Bessel functions: 

Recall the occurrence of modified Bessel functions in §5 Proposition. They are given for 
any complex number \x by the following series: 



for any complex number z in C\R<o- In this section preliminary results on integral 
representations of these functions are reviewed. 

Hankel contours and the Hankel representation of the gamma function: Recall 
that the gamma function at any complex number s with Re (s) > is given by: 



In this section, following [D, pp.225f], the Hankel form of the gamma function is described, 
an integral representation of the reciprocal of T(s) valid for any complex number s. 

First notice that the Laplace transform at any z with Re (z) > of the function x s ~ x on 
the positive real line is z~ s ■ T(s). The following classical result of Laplace now gives the 
inverse of this Laplace transform. For any complex number s with Re (s) > one has: 



for any x > and any fixed £o > 0. Herein deform the path of integration as follows. Fix 
any positive numbers P bigger than R and with R bigger or equal a. Starting at £o — iP 
move parallel to the imaginary axis until £o + iP- From that point move parallel to the 
real axis to the point iP on the imaginary axis. Now follow counterclockwise the circle 
with radius P around the origin. When the parallel to the axis through ia is hit follow 
this parallel until the circle with radius R around the origin is hit. Continue clockwise on 
this circle until R on the positive real axis. Thereafter move to £ — iP m such a fashion 
that the path traced out in the complex plane is symmetric with respect to the real axis. 






dx . 
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The integral of exp (x^)^~ s over the closed path of integration just constructed is zero 
using the Cauchy Theorem. The integrand is rapidly decreasing in the radius |£| of £. 
Letting P go to infinity the above inverse Laplace transform so equals the integral over the 
Hankel contour C a ^: the boundary of the pan in the complex plane of radius R around 
zero with handle of diameter 2a stretching to minus infinity along the negative real line, 



< 


> N. 

ia \ 






-R 

> V 




t >■ 

' R 



Figure A The Hankel contour C a ^R 



on which one comes in from — oo on the branch in {z|Re (z),lm(z) < 0}, passes counter- 
clockwise around zero, and leaves on the branch in {z|Re (z) < 0, Im (z) > 0}. The Hankel 
contour Cr is defined as the limit of a going to zero of the Hankel contours C a ^. 

Returning to the gamma function, the main point now is that on any Hankel contour C the 
absolute value of exp (x£) is exponentially decreasing for any positive x with the absolute 
value of £ going to infinity. Hence one obtains the following Hankel formula: 

valid now for any complex number s, and any positive real number x. 



A Hankel type integral representation of Bessel functions: In this section, fol- 
lowing [WW, 17-231, p. 362], the following Hankel-type integral representations on the 
right half-plane of modified Bessel functions are recalled. 



Lemma: For any modified Bessel function 7 M one has: 

2m JlogC 

for any complex number z with positive real part and any Hankel contour C . 
The key step in the proof is the following integral representation: 

valid for any complex number z in C \ R<o- Indeed, substitute for the reciprocal gamma 
values in the series of I^z) the respective Hankel formulas with x = 1. Interchange the or- 
der of summation and integration. The resulting factor ^ m {z 2 j A.Q'™ 1 j m\ is the exponential 
function at z 2 /4£, as was to be shown. 
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Change variables £ = zr]/2 in the above integral representation of 1^ to get: 

2m J c 

Herein change variables r\ = exp (£ ) to complete the proof of the Lemma. 

The contour log Cr: In the sequel, it is the Hankel contours Cr that are used. For 
a description of the logarithmicalized contour logC^, recall the principal branch of the 
logarithm on the complex plane with the non-positive real axis deleted: 

log(z) = log \z\ + % • arg(z) where — tv < arg(z) < tt. 

Returning to the Hankel contour Cr, the argument of the complex numbers in the upper 
branch of the panhandle in Cr is +7r, of those in the lower branch it is —tv. The elements 
on the circle part of Cr have the form R ■ exp (i9) with — tv < 6 < tv. The contour logC^ 
thus has the following shape: Coming in from plus infinity, move on the parallel through 



log R+in 





1 s» 




l s- 

1 < 



log R—iir 

Figure B The contour logC^ 



the point — in to the real axis to the point logR — in. From this point move up to the 
point log-R -)- iiv on a parallel to the imaginary axis. Finally exit from logR + itv to plus 
infinity on the parallel through itv to the real axis. Notice logC^ C {z|Re(z) > 0} if and 
only if the radius of the circle in Cr is bigger than 1. 

Odds and ends: To conclude, explicit consequences of the above developement are in- 
dicated. First of all using the explicit coordinates discussed above in the logarithmicalized 
Hankel contour with radius R equal to one, the integral representation for 1^ of the Lemma 
specializes to the following Schlafli integral representation: 

I ( z ) = \ I e zcoae cos(n6)d6- [ e - zcosh ^-^ x dx, 

K Jo 77 Jo 

for any complex number z with positive real part [Wa, 6-22]. Direct calculations us- 
ing the series of the respective modified Bessel functions or their Hankel type integral 
representations of the Lemma, prove the recursion rule: 

z ■ I^z) = 2{n + 1) • + z ■ I^ +2 {z) , 

valid for any complex numbers z with | arg(z)| < tv and fx [Wa, 3-71]. 
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The picture is finally completed with Weber's integral: 

/ e~ ax2 x fI+1 lJx)dx= — — — - -e 35 , 
Jo M (2a)^ +1 



valid for any positive number a and any complex number [i with real part bigger than 
minus one. For its proof following [Wa, 13 • 3] develop the Bessel function factor of the 
integrand into its series. Interchange the order of integration and summation. Changing 
variables y = ax 2 , any n-th term of the resulting series is (2a) - ^ +1 ) times the quotient of 
(4a) _n over n factorial times the quotient of f °° exp (— y)y^ +n dy over T(/i+n + l). With 
the real part of \i bigger than minus one, the numerator integral of the third quotient is 
the gamma function at /i+n+1 and thus cancels with the denominator. Thus the series is 
(2a) times the series of the exponential function at the reciprocal of 4a. Applying 
Lebesgue Dominated Convergence then completes the calculation of Weber's integral. 



7. First steps of the Laplace inversion: From paragraph five the Geman-Yor func- 
tions Fqy are for any positive real number a recalled to be given by: 

F GY (a,z) = 



t- (z-2(i/+l)) ' 
where 

D v (a,z)=—— e ^■x» +s -I v ^ T ^{-) dx, 
a Jq \a/ 

for any complex number z with real part bigger than 2(^ + 1). In the sequel the problem 
of inverting this Laplace transform is reduced to the following: 

Lemma: For v bigger than minus one, the inverse Laplace transform at any positive real 
number h of any function Fqy(o, z) is given by: 



where 



u+2 

-hC (2a)^ 
ci = e 1 , 



(z/+l) • e 2a 

where Cr is any Hankel contour with R > Rq, and with integrands absolutely integrable. 



The proof of the Lemma is based on the complex inversion formula for the Laplace trans- 
form reviewed in paragraph zero. 



To prove the Lemma, let z' be any positive real number such that the line {z|Re (z) = z' } 
is contained in the half-plane where Fqy{o,z) is a holomorphic function. Suppose for a 
moment proved that this function is rapidly decreasing, and apply the complex inversion 
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formula to it. Writing out D l/ (a,z), the inverse Laplace transform of Fgy(cl,z) is the 
function on the positive real line given for any positive real number h by: 



e ^ 1 f e h * - X !^ U dxdz. 



W^-ioo " Jo z(z -2(u+l)) 

For an equivalent expression, change variables 77 = 2z+z/ 2 , put zo = 2z' +v 2 , and substitute 
the Hankel-type integral representation of §6 Lemma for the modified Bessel function. 
With c[ = &{v + 1) • (2a) - ^ +4 )/ 2 • ci, the following integral is then to be computed: 

c\ ■ ( tt-7 | / f f e 2a x !/+3 e a cosh ^- ) . e 2 2 -—-dwdxdz. 

1 \2mJ J zo _ ioo J J logCR (z - v 2 )(z - (i/+2)2) 

We claim that for i? sufficiently big the absolute value of the integrand of this integral is 
exponentially decreasing to zero with the absolute values of z, x, or w going to infinity. 

Granting this result, the above triple integral then gives the desired Laplace inverse, and, 
using Fubini's theorem, the order of its integrals can be interchanged. Take the integral 
for Laplace inversion, i.e., the integral over the line {z|Re(,z) = zq\, as inner integral. 
Change variables x = (2a) 1 / 2 t. The Lemma follows on decomposing the denominator of 
the integrand in partial fractions. 

One is thus reduced to prove the last claim about the asymptotic behaviour. For the 
calculations recall | exp (£)| = exp (Re (£)), for any complex number £. For the asymptotic 
behaviour in the absolute value of w, reduce to elements w = x ± iw in log Cr. For these 
cosh(w) = — cosh(x). Hence the absolute value of the hyperbolic cosine factor of the nu- 
merator equals exp (—(x/a) cosh(Rety)) whence the required asymptotic behaviour in the 
absolute value of w. 

For the asymptotic behaviour in z, notice | exp (— wz 1 / 2 )] = exp (— Re (wz 1 / 2 )) and recall 
z 1 ! 2 = exp ((1/2) • (log \z\ + i arg(z))). The argument of z converges to 7r/2 with \z\ going 
to infinity. Now Re (wz 1 / 2 ) equals \z\ 1 ^ 2 times Re (w) cos(arg(z)/2) — Im (w) sin(arg(z)/2). 
Herein the cosine is positive and bigger than cos(7r/4). Thus Re (wz 1 / 2 ) is positive if the 
real part of any w is big enough. Hence choose R big enough to have the desired asymptotic 
behaviour in the absolute value of z. 

The asymptotic behaviour in x is determined by exp (— x 2 /2a), thus completing the proof. 

8. Computation of certain Laplace transforms: Lacking a suitable reference, this 
section computes the inverse Laplace transforms identified in §7 Lemma. 

For any a and (3 in C, consider the functions on the positive real line given by: 

9a At) =^ (e +a P Erfc + pVT^J + Erfc - pVt^j ) , 
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for any positive real number t. Then one has the following two results: 
Lemma: If the real parts of a and a 2 are positive, one has: 

L (f a ,p) (Z) 



g-a\/z" 



for any complex number z in C \ R<o with real part bigger than |Re 
Corollary: If the real parts of a and a 2 are positive, one has: 

L (g a ,p) (z) = - _ /32 , 

for any complex number z in C \ R<o with real part bigger than Re ((3 2 ). 

The Corollary follows from the Lemma upon decomposing the denominator in partial 
fractions and using the linearity of the Laplace transform. 

We use the following two results proved mutatis mutandis in [D, Beispiel 8, p.50f]: 
Lty«)(z) = e~ aVT where ^(t) = -2= • H • e _fi , 



L(xa)(z) = — 7=- where x<*(t) = • t~ * • e ! ? 



for any complex number z in C \ R<o and any positive real number t. 
Subtracting \ a from f a ,/3, the proof of the Lemma reduces to show the identity: 



L (-(3 ■ e ° W 2 -* Erfc (fly/t + -^=j ) (z) = -0 



Multiplying any nice function / with exp (a-) induces a shift by —a in its Laplace trans- 
form: L (exp (at)f(t))(z) = L (f)(z — a). Using this with a = (3 2 one is further reduced to 
calculating the Laplace transform of the above complementary error function factor only. 

Since Re (a) is positive, the real part of (3v}l 2 + (a/2)u~ 1 / 2 goes to plus infinity with u 
converging from the right to zero. The Fundamental Theorem of Calculus thus gives: 

/ \ 2 

t 
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Using the transform-of-an-integral property L (J** f(u) du){w) = w 1 L (f)(w) the Lapla- 
ce transform of this complementary error function at w = z — (3 2 is given by: 



2(3 



+ 



2 

2(z-p 2 y 



a 



1 -U^-s^) 

. t -2 . e v 2 vw 



t 2 • e 



(z-/3 2 ). 



Using the Laplace transform of Xa, the first Laplace transform of this sum equals 



Using the Laplace transform of ip a , the second Laplace transform of this sum equals 

The above sum hence equals exp {—af3) exp (— a^/z~) / (z 1 / 2 (z 1 ^ 2 + /?))• The identity to be 
proved follows upon substituting this last expression. This completes the proof. 

9. Two intermediate results: In this section the Laplace inversion of §7 is resumed 
concentrating on the two summands of the expression derived in §7 Lemma. 

If v> — 1, for any real numbers a, h > and 6, thus define more generally: 



i 

2~kI 



logCR JO 



e ~ x2+x ^ cosh (™)^+3 L -i 



-W\fz 



h 



dx dw 



with log(7R the logarithmicalized Hankel contour discussed in §6. For these integrals one 
has the following two results, the first of a more structural and the second of a more 
technical nature: 



Lemma: If u>— 1, for any real numbers R>1, and a, h > 0, and b one has: 

cosh(w) 

{ '-'h\'i) r r b [ji) ){tc) • n_ f , H h [ - 

'log C R 

where on both sides the respective integrands are absolutely integrable. 



1 l{h) = C2-T{v + A)~ i f^ ^ {F b (h) + F_ b (h))(w).H_ { ^ ) [- C ^^l)dw, 



Corollary: // u> — 1, for any real numbers p = \ogR > 0, and a, h > 0, and b one has: 



with Re (£) and Im (£) the real respectively the imaginary part of any complex number £. 
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To explain the notation in these two results, C2 denotes the following constant: 

The Hermite functions of degree \i are discussed in paragraph zero, and the two other 
functions that occur are given by: 

F b (h)(w) =e- b -Erfc (-^= + ^V^j , 

^(/ i )H=e X ^ COShM (F & (/,)H + F_ b (/ l )H), 
for any complex number w. 

Proof of the Lemma: For the proof of the Lemma choose a Hankel contour Cr with 
R so big that for any element w in log Cr also the real parts of w and w 2 are both posi- 
tive. At any w on logC^ substitute in I%(h) the inverse Laplace transforms calculated in 
§8 Corollary. Interchange the order of integration using the absolute integrability of the 
integrand. This gives the expression of the Lemma for l%(h). The integrand herein is a 
holomorphic function on C \ R<o- Using the Cauchy Theorem, the value of the integral 
is independent of the Hankel contour Cr chosen as long as logR is non-negative. This 
completes the proof. 

Proof of the Corollary: Given the Lemma, the proof of the Corollary is an exercise in 
path integration. Put G x = G a x b (h). Changing the order of integration in the Lemma, 



/ e-*V +3 / G x (w)dwdx 

Jo 2lvt JlogC R 



f oo 

C2 



' log Cr 

is to be calculated. Concentrate on its inner integral and consider the following subpath 

P = P C + P 00 with P c = P Ci+ - P c ,_ 

of the path log Cr. Here the path P c ,- starts from logR and moves parallel to the imagi- 
nary axis to the point logR — m. The path P c ,+ starts from logR and moves parallel to 
the imaginary axis to the point logR + in. The path P^ moves from logi? + iir parallel 
to the real axis to +00. For any x > 0, the inner integral thus breaks up as follows: 

h L ( GAw) ~ GAw) ) dw+ inI P GAw) dw - h j P GAw) dw ' 

with w the complex conjugate of w. Using the series expansion of the exponential and 
the complementary error functions, G x is compatible with complex conjugation, i.e., G x 
evaluated at the complex conjugate of any complex number w is the complex conjugate of 
G x at w: 

G x (w) = G x (w). 
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Write the elements of P^ as w = y + in with y > log R, notice dw = dy and change vari- 
ables. Using the compatibility of G x with complex conjugation, one obtains the improper 
integral from log R to infinity of the imaginary parts of G x (y + in) . 

On the circle part P c of P one has to be a bit more careful about the volume forms. The 
elements of P C) + are parametrized by w = \ogR + id with < 6 < n, whereas those of 
P c - are parametrized by logR — i9 with < 9 < n. Changing variables accordingly and 
integrating from zero to 7r, the induced volume form in the P C; + -integral is dw = i d6, 
whereas that on the Pc^-part is dw = —i d0. Using the compatibility of G x with complex 
conjugation, this completes the proof of the Corollary. 

10. The final explicit calculations: The aim of this section is to explicitly compute 
the integrals of §9 Corollary thus essentially proving the valuation formula §4 Theorem. 

Recall the integral l%(h) of paragraph nine, the functions Erfc , introduced in para- 
graph zero, and the function E b (h) introduced in paragraph four. In c^ih, b) of paragraph 
nine divide by n, add a gamma factor, and drop the exponential in h and b to obtain: 

c 3 = c 3 (v) = ^ -r(z/ + 4). 

In view of §7 Lemma, the following result is the key for proving the valuation formula: 

Lemma: If u>— 1, for any real numbers a, h>0 and b, the function l^(h) is given by: 

2 • c 3 • e h * ■ r H, v+A) (-^B-) cos(^) d9 
Jo V v2a / 

+C8 • A . e 3>£ . £ #_ ( „ +4) (^^) (ey b E b (h)(y) + e-y b E_ b (h)(y))dy. 

The Lemma is proved by computing the two integrals of G x :— G^ b of §9 Corollary upon 
choosing R = 1 there. Start with the following integral: 

poo 

/ Im (G x )(y + in) dy . 
Jo 

Abbreviate j3± = y/V2h± (b/2)\/2h and notice cosh(y + in) = — cosh(y). Write the 
complementary error functions occuring in G x (y + in) as improper integrals starting from 
zero as in paragraph zero. Multiply out the expressions of the exponents to obtain: 

e Mv+t*)b ET f c (p +i _? = \ = eTi . e ±vb e -(u+p±) . e V h Jdu. 

V V2h J Jo 

Since j3± is a real number and u can be taken as real numbers, the imaginary parts in G x 
are determined by the imaginary parts of the exponentials in the integral. These are given 
by the functions E± b (h) as indicated, completing the calculation. 
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As a next step, 

Re (G x )(iO) dO 



o 



is calculated. In contrast to the above argument, in the complementary error func- 
tions occuring in G x (i0) now the following paths of integration are used: Abbreviating 
(3± = ±(6/2)v / 2/T, first move from (3± + i9 to /3±, then continue from (3± to plus infinity 
along the real line. Since cosh(i6>) = cos(6>), the above integral then equals: 

£ e ^-cos(e) | ( Erfc + Erfc cos{%) + J^^ + + m, 

upon abbreviating for any angle 9: 

P±+- " 

P± 



i±(0) = -Re ( e ±l6b / V2h e"^ du 



In the first of these two last integrals notice Erfc (/?+) + Erfc (/?_) = 2 since /3_ is minus 
(3 + . To calculate <f>± change variables u = (3± + iw/V2h in the integral. Write the factor 
ij\p2~h that is picked up as exp /2)) / 'y/2h '. Multiplying out the expression obtained in 
the exponent, it follows 



Jb± y2h Jo 



e 

>(3± 



The real part of this expression is determined by the real part of the exponential functions 
in the integral. Abbreviating x u = b(9 — u), the values of the cosine at 7t/2±x u thus appear 
as factors. A shift by tv/2 turns a cosine into a sine as follows: cos(-7r/2 ± x u ) = =F sin(x u ). 
Thus 4>- is minus 0+, completing the proof. 



11. First part of the proof of the valuation formula: The proof of the valuation 
formula of §4 Theorem is in two steps. As a first step, §4 Theorem is in this paragraph 
established for v bigger than minus one. As a second step, this equality is extended in the 
next paragraph to any complex number v using analytic continuation. 

Thus let v be bigger than one. Recalling §7 Lemma, the valuation formula of §4 Theorem 
then is obtained by subtracting 11(h) from l q v+2 (h) and thereafter multiplying this differ- 
ence with the constant c\. Substitute the expressions computed in §10 Lemma. Using 
exp {h{y + 2) 2 /2) = exp (hu 2 /2) exp (2h(u + 1)) with the trigonometric summands one is 
reduced to show that 

C t * rig = £ H_ {v+A) (-^) cos(irf) d6 

upon multiplication with c gives the first two terms in the sum of §4 Theorem. As a 
first step, write C t * rig as a double integral using the defining integral representation of the 
Hermite function factor of its integrand. Applying Fubini's theorem and interchanging the 
order of integration the integrand of its inner integral is exp (2x(2q)~ 1 / 2 cos 9) ■ cos(z/6>). 
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Substitute for this last integral using Schlafli's integral representation of paragraph six. 
Then reverse the order of integration to get: 



7 + sin(z/7r) 



poo 

Jo 



( cosh(x^ 



where 



/•OO 

J = / e"*' 
The proof thus reduces to show: 



2a; 



I=--e^-{2q) ~^~(l + 2g> + l)). 
Indeed, abbreviating a = q/2, change variables w = a~ x l 2 x in / to obtain: 



a 2 + 2 / 
Jo 



Apply the the recursion relation for modified Bessel functions of paragraph six 
2(i/ + 1) 

Jo 



I = a 



POO POO 

2(i/ + 1) / e- aw2 w u+2 I I/+1 (w)dw+ / e" au,2 «;" +3 ^ +2 («;)rf«; 



The evaluation of / is thus reduced to evaluating two Weber's integrals as reviewed in para- 
graph six. Substituting their respective values, the above identity follows. This completes 
the first step of the proof of §4 Theorem 

12. Second part of the proof of the valuation formula: This second part of the 
proof of §4 Theorem extends its validity from v bigger than minus one, as established in 
the previous paragraph, to v any complex number. 

This reduces to show the following two results. First, the normalized price is an entire 
function in v. Second, the right hand side of §4 Theorem is a meromorphic function on the 
complex plane. Indeed, these two functions agree on real numbers v bigger than minus one. 
Using the identity theorem they so agree on the complex plane as meromorphic functions. 
With one of them entire, the other one is entire, too. 

In Yor's notation A^\h) = f£ exp {{2{B U + vu)) du recall CM(h,q) = E[f(AM(h))], 
where / is given by f(x) = (x — q) + , for any real number x. The proof of this normalized 
price being entire in v is then further reduced to show 



E 



f(A h ) ■ e 



an entire function in v where A is the process A(°). Indeed, this follows using the Girsanov 
identity E[f (AM (h))] = exp (-z/ 2 h/2) ■ E[f(A(h)) ■ exp(uB(u))] of [Y, (l.c), p.510]. 
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For proving the above function entire develop the factor exp(uBh) of f(Ah) ■ exp(i/Bh) 
into its exponential series. Suppose computing the expectation of the resulting series term 
by term is justified for any complex number v. For any complex number is, this then gives 
a power series in v and thus explicitly shows E[f(Af l ) ■ exp(uBh)] holomorphic at v. 

Interchanging the order of integration and summation is justified using Lebesgue Domi- 
nated Convergence if the following is true. The expectations of the absolute value of any 
single term of the series for j{Ah) ■ exp [yBh) exist and the series so obtained converges. 
Using the Cauchy-Schwarz inequality this is implied by the series: 



1/2 



being convergent for any complex number v. Implicit herein is that j{Ah) is square inte- 
grable. This is implied by Ah being square integrable. Any n-th moment in particular of 
Ah has been computed in [Y, (4.d"), p. 519] as: 



(~l) n 2 y^ (-l) n - k -h- 
(n!) 2 ^k=o (n-k)\ ■ (n+k)l 



In particular the second moment of Ah is thus finite, as was to be shown. What regards 
the even order moments of the Brownian motion B at time h, they are computed as: 

E[Bl n ] = W. r( „ +1/2)i 

for any non-negative integer n. The above series thus converges using the ratio test. This 
completes the proof of the normalized price being an entire function in the paramater v. 

The proof of the five-term sum of §4 Theorem being a meromorphic function in v is based 
on the Hermite functions being entire also in their degree. Indeed, from [L, (10.2.8), p. 285] 
one has for any complex numbers [i and z the representation: 

2*-r(l/2) ^1 Ai 1 ,\ 2"-r(-l/2) ^l-f< 3 , 



H - {z)= T( ( l-,)l2) ^\--2-2' Z ) + r(-„/2) 

Herein the reciprocal of the gamma function is entire by construction, and the confluent 
hypergeometric function $ is entire in its first and third variable. 

There is a localization principle for proving entire a function on the complex plane. Indeed, 
one has to show analyticity at any fixed complex number. For this, one can restrict the 
function to any relatively compact or compact neighborhood of this fixed complex number. 

Apply this localization principle to the factor of the trigonometric term: 



and restrict v to belong to any suffiently small compact neighborhood U of any point fixed 
in the complex plane. The above integrand is analytic as a function in v and smooth as 
a function in 9. Thus the absolute values of its derivatives with respect to v of any order 
are bounded on the product of U and the closed interval between zero and ir. In this local 
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situation, using a standard consequence of Lebesgue Dominated Convergence, differenta- 
tion of A with respect to v is thus by partial differentiation with respect to the parameter 
v under the integral sign. This proves A entire as a function in v. The trigonometric term 
Ctrig is obtained by multiplying A with c. Dividing off from c the gamma factor T(z/+4), 
it has a removable singularity in v = —1. Thus C tr i g is meromorphic in v with at most 
simple poles in the integers less than or equal to minus four. 

What regards the hyperbolic terms, the claim is that as a function of v they can be ex- 
tended as analytic functions from v bigger than minus one to to the whole complex plane. 
For this question it is sufficient to consider the function B given by: 



for any real number y. Extension of this function is essentially by reduction to the case 
where the real part of the degree — (z/+4) of the Hermite function factor in the integrand 
of B{y) is negative, or equivalently, the real part of v is bigger than minus four. 

To fix ideas first consider the case where v is such that the degree of the Hermite function 
factor in the integrand of B(u) is a non-negative integer. This Hermite function is then 
the corresponding Hermite polynomial. In the integrand of B(u) the absolute values of 
the Hermite function factor and the exponential function factor so have linear exponential 
order in the variable y. The decay to zero of the absolute value of the function E b ^, how- 
ever, is of square exponential order in y. It thus dominates the asymptotic behaviour with 
y to infinity. The absolute value of the integrand of B(u) is so majorized by an integrable 
function, and B can be extended to the above values of v. 

For the general case of the reduction, fix any v of real part less than or equal to minus 
one that is not an integer. Apply the above localization principle and let v belong to a 
sufficiently small compact neigborhood U in the half-plane {Re (z) < —1}. Shrinking U if 
necessary assume that for any v in U the degree of the Hermite function in the integrand 
of B(u) is not an integer. Using the recursion rule for Hermite functions of paragraph 
zero express the Hermite function factor of B(u) in terms of weighted Hermite functions 
of negative degrees. Further shrinking U if necessary, assume that the so obtained relation 
represents B on U. Herein the weighting factor for the respective Hermite functions are 
given as powers of z = (2q)~ 1 ^ 2 cosh(y) times polynomials in v. The absolute values of 
the polynomials in v can be majoriozed uniformly on U. The problem thus reduces to 
majorize by an integrable function on the positive real line in the variable y finitely many 
functions on U times the positive real line sending v and y to: 



where k, I range over finitely many non-negative integers and k is such that u + 4 + k is 
positive. The leading terms of the asymptotic expansion of paragraph zero for any Hermite 




where b(v) equals ±u or ±(z/+4) and with: 





Valuation of arithmetic-average options 



25 



function H^(z) with degree \i any complex number with negative real part has order . 
Asymptotically with y to infinity, the Hermite function with the smallest positive number 
v+4+k thus dominates the other Hermite function factors in the above functions. With v 
ranging over a compact set, there is a minimal such degree on U. Similarly, there are such 
majorizing choices £* for the factors cosh^(y), and v* for the absolute values of the factors 
exp (b(u)y) , and u** for the absolute values of the factors E b ^(h)(y). A four-factor- 
majorizing function on the positive real line thus results whose asymptotic behaviour with 
y to infinity is governed by the square-exponential decay to zero of the corresponding 
factor E b ^**^(h) and which is integrable. 

If the real part of v is bigger than minus one, the above argument holds in a simplified 
form. The upshot so is that any complex number not an integer has a sufficiently small 
compact neighborhood such that the absolute value of the integrand of B{y) on U times 
the positive real line can be majorized by an integrable function on the positive real line. 
Herein, compact neigborhoods can be replaced by relatively compact neigborhoods mu- 
tatis mutandis. Using Lebesgue Dominated Convergence, B can thus be extended as a 
continuous function to the whole complex plane with the integers less than or equal to 
minus four deleted. 

The idea for showing B analytic as a function of v on the complex plane with the integers 
less than or equal to minus four deleted is as follows. Show that differentiation of B is by 
differentation under the integral sign and use that its integrand is entire as function of v. 
For this again first localize to v in any sufficiently small compact neighborhood containing 
no integers less than or equal to minus four. The aim is then to majorize the absolute value 
of the derivative with respect to v of the integrand of B by an integrable function indepen- 
dent of v as above. The above argument for getting such a majorizing function is based 
on a comparison of decay rates. The integrand of B has one factor which on the positive 
real line decays to zero of square exponential order whereas the other factors explode of at 
most linear exponential order. This situation is preserved on differentiation with respect 
to the parameter v. In particular, differentiating with respect to the degree the asymptotic 
expansion for Hermite functions on the right half-plane gives an asymptotic expansion for 
this function's partial derivative with respect to the degree. 

At this stage it remains to extend B analytically to the intgers less than or equal to minus 
four. However, B remains bounded in any punctured compact neigbourhood of such an 
integer. Thus B can be extended to an entire function, completing the proof of §4 Theorem. 



Part III 

13. Remarks about hedging: This paragraph's aim is to compute the Asian option's 
Delta and discuss how the seller's hedging portofolio is determined by it. 

Delta, and similarly the other local hedging parameters, are computed by partially differ- 
entiating the price function Ct,T(K) of paragraph three: 

C t , T (K) = e-^) . St ■ a2{T 4 _ to) ■ C^(K q). 
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To simplify notation, in the sequel as many arguments of a function as possible are sup- 
pressed, thus writing C t = Ct,r(K) and C^ v > = C^ u \h, q) in particular. 

Hedging of the Asian option as a particular case of the general theory of hedging European- 
style contingent claims in a complete Black-Scholes economy has been discussed in [K, 
p.23f]. The seller's hedging portofolio IT is determined using the martingale representation 
of the conditional expectation of the discounted pay-out of the Asian option at its time of 
maturity, i.e., using the following stochastic differential equation: 

e -r{t-t ) Ct T{K) = c kuT (K) + a f e-^ s -^Ii s dB s . 

Jt 

Apply to the martingale of the left hand side the ltd formula. Comparing diffusion coeffi- 
cients, the hedging portofolio Il t at any time t between t and T thus is given by: 

dC 

U t = Sf A* where A, 



dS t 

The option's Delta A t is computed in the case where q is positive and §4 Theorem applies. 
The partial derivative of q with respect to St times St being minus q, one has: 

dS t a 2 (T-t )V [n ' q) Q dq [n,Q) )- 

One is thus reduced to computing the partial derivative with respect to q of C(") at (h, q). 
Think of this as the option's normalized Delta. Further reduce as follows. Write the 
five-term sum of §4 Theorem for CM in the form: 

CM = d ■ Arig + d ■ D hyp , 

setting: 

, ,, , I>+4) (2g)^ 
and where the modified trigonometric term -Dtrig and the modified hyperbolic term: 

e • -Dhyp = Chyp,^ + 2 + C hy p ) _( I , + 2 ) - Chyp,^ - Chyp,-^ 

are functions in the variables v, h and in q = q(h). The partial derivative of C<M with 
respect to q then is: 

Thus, one is further reduced to compute the partial derivatives with respect to q of the 
above functions d, -Dtrig, and Ch yp ,b with b equal to ±v or ±(z/ + 2). The first of these is: 
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To compute the partial derivatives of the other two functions recall from paragraph zero 
the recursion formula for the derivative of Hermite functions. Granting for a moment 
that one can justify the differentiation under the integral sign, the partial derivative with 
respect to q of .Dtrig at (v, h, q) is: 

With the functions Ej,(h) of paragraph four independent of q, the partial derivative with 
respect to the variable q of any Ch yp ,6 at (v, h, q) similarly is: 

2 £ 2(z/+4) f°° TT /cosh(y)\ yb „ , IWXI 

To justify the differentiations under the integral sign first notice that differentation is a 
local concept and thus the parameters v, h can be restricted to vary in a fixed compact 
set not containing points with q equal to zero. The integrands of -Dtrig and -Dh yp then are 
integrable and differentiable functions not only in the variable 9 respectively y but also 
in the variables v, h, and in q. Consider the maxima in h of the absolute value of the 
partial derivatives with respect to q of -Dtrig and -Dh yp , for any fixed triple of the other 
variables. In particular with q bounded and bounded away from zero the so obtained 
functions in the variables 9 respectively y, v then are majorized by integrable functions. 
A standard application of Lebesgue Dominated Convergence thus completes the argument 
for the differentiations under the integral sign, and thus completes the calculations. 
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